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Abstract. We investigate the design and functionality of a network of loop-shaped 
charge qubits with switchable nearest-neighbor coupling. The qubit coupling is 
achieved by placing large Josephson junctions at the intersections of the qubit loops 
and selectively applying bias currents. The network is scalable and makes it possible to 
perform a universal set of quantum gates. The coupling scheme allows gate operation 
at the charge degeneracy point of each qubit, and also applies to charge-phase qubits. 
Additional Josephson junctions included in the qubit loops for qubit readout can also 
be employed for qubit coupling. 



PACS numbers: 74.81.Fa, 03.67.Lx, 85.25.Hv, 85.25. Cp 



Superconducting qubit network with controllable nearest neighbor coupling 



2 



1. Introduction 

During the last six years it has been experimentally proven that superconducting circuits 
can serve as quantum mechanical two-level systems, qubits, to be used for quantum 
information processing [H El El El El E] • Besides the experiments with individual qubits, 
several experiments have been performed so far on two permanently coupled qubits 
[3 El El HOI \U\ For instance, to observe the coupling of two charge qubits, the qubit 
islands have been permanently coupled via a capacitor, and the strength of the coupling 
has been varied by tuning the qubits in and out of resonance with each other (by varying 
the gate voltage) jZUHj- 

In order to build a functional, scalable quantum computer, a network design is 
needed that allows coupling of an arbitrarily large number of qubits, with the possibility 
to switch on and off the coupling by means of external control knobs. In principle, 
coupling of only nearest neighbor qubits is sufficient to perform a universal set of gates 

PI- 

Theoretical schemes for variable coupling of charge qubits have been intensely 
discussed in literature. Couplings via inductive and capacitive elements have been 
examined as well couplings via linear LC-oscillators and Josephson junctions ITU 
[T71 HE] . A standard approach to achieve a variable coupling is to employ a SQUID- 
type geometry either for the qubits |18 J or for the coupling element [T5] . to be able 
to control the Josephson energy by an external magnetic flux. A somewhat different 
approach has been suggested in [20], where the qubits are coupled via another charge 
qubit thus creating a variable capacitive coupling. 

Recently, a different way to achieve a variable inductive coupling has been 
suggested, namely to let the charge-qubit loops intersect and share a coupling Josephson 
junction or SQUID. The interaction is then controlled either by varying the magnetic 
flux in the qubit loops (or the coupling SQUID) [21], or by applying bias currents to 
the coupling Josephson junction 

In this paper we give a detailed analysis of the qubit network based on the coupling 
method proposed in [22] • The idea of this method is to couple loop-shaped charge qubits 
by letting the circulating loop currents, which are sensitive to the charge state of the 
qubit island, interact. This is done by placing a non-linear oscillator - a large Josephson 
junction - at the intersection of the qubit loops. Such a coupling can be made variable by 
using the fact that in the absence of an external magnetic field, the persistent currents 
in the qubit loops are absent (for symmetric qubits with equal Josephson junctions). 
However, when a dc current bias is applied to a coupling Josephson junction the 
symmetry is broken and currents start to circulate, the magnitude of the currents being 
dependent on the bias current. These currents interact with the oscillator, resulting in 
a variable effective qubit-qubit coupling. A similar coupling effect can be accomplished 
by inserting large readout Josephson junctions in the qubit loops [2 . As we show in 
this paper, applying current to one of the readout junctions allows one to measure the 
state of the corresponding qubit without disturbing the other qubits; however, when 
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two neighboring readout junctions are biased, the qubit-qubit coupling is switched on. 

The advantage of the current-biased coupling scheme is that it does not need any 
local magnetic fields to control the coupling, fields which could create unwanted parasitic 
long range interactions. An important feature is the possibility to operate at the qubit 
charge degeneracy point, where the decoherence effect is minimized and where 

the gate operations are very simple. This coupling scheme can be also extended to 
charge-phase qubits which are still less sensitive against decoherence due to flatter band 
structure [2j. With this coupling mechanism, neighboring qubits in an arbitrarily long 
qubit chain can be coupled, and several independent two-qubit gates can be performed 
simultaneously. The fundamental entangling two-qubit gate is a control-phase gate, 
which together with single-qubit gates constitute a universal set of operations. 

The structure of the paper is as follows: In section |2] we explain the principles of 
the coupling by considering the simplest case of two coupled qubits. We estimate the 
maximum coupling strength, evaluate the residual parasitic couplings, and investigate 
the charge-phase regime for the qubits. In section 01 we add measurement junctions to 
the qubit circuits and investigate how to use these junctions as read-out devices and the 
means to create qubit coupling. In section 0] we generalize the derivation of section El to 
a multi-qubit network with an arbitrary number of coupled qubits. Finally, in section 
EJ we discuss how to use the network for quantum computing. 

2. Controllable coupling of two qubits 



Qubit 1 



Qubit 2 



—I— r —I— r 



J V gi 



I 
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Figure 1. A system of two coupled charge qubits. V g i control the individual qubits 
whereas lb controls the coupling of the two qubits. 



To more clearly explain the principle of the qubit coupling, we first consider the case 
of two coupled qubits. The qubits consist of single Cooper pair boxes (SCB) with loop- 
shaped electrodes. 01211 • To create coupling between the qubits, a large-capacitance 
Josephson junction (JJ) is placed at the intersection of the qubit loops, see Fig. [U The 
physics of the coupling is the following: As long as no magnetic flux is applied to the 
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qubit loops, and no bias current is sent through the coupling junction, there are no 
circulating currents in the qubit loops. However, when the bias current is switched 
on, circulating currents start to flow in the clockwise or counterclockwise direction 
depending on the charge state of the Cooper pair box. These currents displace the 
coupling junction oscillator and change its ground state energy, leading to an effective 
qubit-qubit interaction. The strength of the interaction is proportional to the bias 
current through the coupling junction. In the idle state, when the bias current is 
switched off, small phase fluctuations at the coupling junction generate permanent 
parasitic qubit-qubit coupling. However, this parasitic coupling can be made small 
compared to the controllable coupling by choosing the plasma frequency of the coupling 
junction uj b = \J 2E h jE h c to be small compared to the Josephson energy Ej. This 
requirement implies that the junction charging energy E h c must be small, 

E b c <^u b ^E b j, (1) 

i.e. the coupling junction must be in the phase regime. This is the most essential 
requirement for the qubit coupling under consideration. 



Ti CJ { Ze \ 

-Li.i = I Vli + 02i) + 2 / 2 ? 2 ( y V 9i ~ + E J ( COS + COS fci) , (3) 



2.1. Circuit Hamiltonian 

We begin the evaluation of qubit coupling with the derivation of a circuit Hamiltonian. 
To this end, we first write down the Lagrangian L of the circuit in Fig. ^ The 
Lagrangian consists of the respective Lagrangians of the SCBs and the coupling JJ, 

2 

L — Lscba + Ljj . (2) 
i=i 

Assuming the single Cooper pair boxes to consist of identical junctions with capacitance 
C and Josephson energy Ej, and following the rules described in e.g. Refs. [221120], we 
write the corresponding Lagrangian on the form, 

where (f>u (<p2i) is the phase difference across the left (right) Josephson junction of the i- 
th SCB, and C g is the gate capacitance. The Lagrangian of the coupling JJ includes the 
electrostatic energy and the Josephson energy of the junction, and also the interaction 
energy of the junction with applied bias current lb, 

where ip is the phase difference across the coupling junction. 

The flux quantization condition in each of the qubit loops allows the elimination 
of one of the qubit variables from the Lagrangian. We assume that there is no external 
magnetic flux in the loops since magnetic flux will not be used to control either the 
qubits or the qubit interaction, and we also assume that the loop self-inductances are 
negligible. Then the flux quantization equation takes the form, 

+ ip = 0, +; 2 -(p = 0. (5) 
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where we introduced new qubit variables, 



4>+,i = 4>u + fa- (6) 



By virtue of relations the gate capacitance terms will take the form 

(h 2 C g /2(2e) 2 )(2eVgi/h + ± 0/2) 2 introducing a capacitive interaction of the SCB 
with the coupling junction. From here on, the upper (lower) sign corresponds to the 
first (second) qubit. Similarly, the appearance of the variable (p in the SCB Josephson 
terms introduces an inductive interaction between the SCB and the coupling junction. 

At this point, we are ready to proceed to the circuit Hamiltonian. By introducing 
the conjugated variables, n; = (l/h)(dL/d(j)- t i), and n = (l/h)(dL/d(p), which have 
the meaning of dimensionless charges (in the units of Cooper pairs) on the SCB 
and on the coupling JJ, respectively, and then applying the Legendre transformation, 
H = J2ih n i4 ) -,i + hrnp — L, we get, 

H = j^H SC B,i + Hjj + H c . (7) 

i=\ 

Here 

H S CB,i = E c (m - n gi f - 2Ej cos - cos (8) 

is the SCB Hamiltonian, where E c = (2e) 2 /2C*s, = 2C + C g is the total capacitance 
of the qubit island and n gi = C g V gi /2e is the (dimensionless) charge induced on the 
qubit island by the gate voltage. The JJ Hamiltonian is 

Hjj = E b c (n- 2 _ e* cos ^ _ l Ih(pi (Q) 

where E b c = {2e) 2 /{2C b + C s ). The last term in Eq. (0), 

He = ^fb (n - %1 ~ ng2 ) ((n 2 - n g2 ) - (m - n gl ))-^E b c ( ni -n gl )(n 2 -n g2 ), (10) 

describes capacitive interaction of the qubits and the JJ, and also direct qubit-qubit 
coupling, induced by the gate capacitance. 

The Hamiltonian (J7J) is quantized by imposing the canonical commutation relations, 
[4>^j,n k ] = iSjk, [<p,n] = i. To incorporate the Coulomb blockade effect, we take 
advantage of the periodic SCB potential and impose periodic boundary conditions on 
the wave function with respect to the phase This results in charge quantization 

on the island. Keeping the system at low temperature {ksT < Ec) and close to the 
charge degeneracy point n g = 1/2, restricts the number of excess charges on the island 
to zero or one Cooper pair. Assuming the charge regime, Eq ^> Ej for the SCB, and no 
transitions to the higher charge states to occur during qubit operation, we truncate the 
SCB Hilbert space to these two lowest charge states. Then the single qubit Hamiltonian 
reads, 

Hi = (1 - 2n gi ) a zi - Ej cos ^a X i, (11) 
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while the capacitive interaction QlUjl takes the form, 

He = ^E h c (n - ^) (a z2 - a zl ) - ^E b c a zl a z2 . (12) 

The qubits and J J interact both capacitively, Eq. (fT2j) . and inductively (the last term 
in Eq. pip). It has been noticed by Shnirman et al. that the capacitive interaction 
can be fully transformed into an inductive one. This can be done by using a unitary 
rotation conveniently combined with a gauge transformation eliminating the gate charge- 
dependent terms in Eqs. (JTSJ) and (this is possible since the charge on the coupling 
J J is not quantized). The corresponding unitary operator is 

n gl -n g2 \ _ C g 



U = exp \-ia{o z2 - a zl )(p] exp \ i tpj , a = (13) 

It is straightforward to check that the transformed part of the Hamiltonian, U'(Hjj + 
Hq)U, does not contain any interactions, while the whole interaction is concentrated in 
the Josephson term of the qubit Hamiltonian, 

U j HiU = (1 - 2n gi ) a zi - Ej cos ^ [cos(2ay?) a xi ± sin(2a^) a yi ] . (14) 
The a-dependent correction is small when the gate capacitance is small, C g C%. 

2.2. Controllable qubit coupling 

Let us consider the main part of the inductive interaction. Making use of assumption, 
Eq. (JTJ), we consider small phase fluctuations across the JJ, 7 = — (po, around the 
minimum point ipo, and expand the (^-dependent terms in Eq. ()14)) in powers of 7 <C 1. 
The quantity <p Q is determined by the applied bias current 

Sm ^ 0= 2^- (15) 
To zeroth order with respect to 7 we get a free qubit Hamiltonian which, after additional 
rotation U' = exp [ia{a z2 — er z i)<^o]) takes the form, 

Hi = ~y (1 - 2n gi ) a zi - Ej cos yff^. (16) 

Controllable qubit-qubit coupling results from the terms which are linear in 7 in the 
expansion of Eq. (|14|) . Neglecting the effect of the small a in these terms (which will 
be considered in the next section), we obtain, 

H int = -Ej sin — j(a x i + a x2 ). (17) 

It is convenient to combine these terms with the quadratic potential of the JJ that 
approximates the tilted Josephson potential near its minimum. We then get the total 
Hamiltonian on the form, 

H = J2 H — XEj - - a x ia x2 + H osc (18) 

i 4 cos (f 
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where A = Ej/Ej, and 

(19) 



1 • ' 2 

if osc = E b c n 2 + -Ej cos y? 



, sin(o9 /2) , , 

7 + A - VV ^ U/ ' ((T x i + £7,2 
2 COS </9 



is the Hamiltonian of a displaced linear oscillator associated with the coupling JJ. 

The further analysis is significantly simplified if one assumes the qubit-oscillator 
interaction to be small, A < 1, and the oscillator to be fast on a time scale of qubit 
evolution, u;& ^> Ej. In fact, these assumptions are not needed when the qubits are 
parked at the charge degeneracy point, n gi = 1/2, because in this case the interaction 
term commutes with the qubit Hamiltonians, and the problem is exactly solvable |27j . 
The assumption on A can be relaxed for the two-qubit circuit, however for a multi-qubit 
network it becomes essential as discussed later. 

The imposed constraints together with Eq. lead to the following chain of 
inequalities: 

Ej < w fe < E b j. (20) 

Under these constraints, one can neglect the excitation of the oscillator, which at low 
temperature will remain in the ground state. This does not significantly differ from 
the ground state of the Hamiltonian in Eq. (|T9"j) . For instance, the estimate for the 
amplitude of the first order correction reads, 

EcQ. - 2n g ) Ej 

C °^ 1 ~ (^)V2( cosv , )3/4 ~ * A ' 

Therefore one can average over the oscillator ground state and drop the oscillator energy 
term, because it does not depend on the qubit state configurations. 

Summarizing our derivation, after integrating out the oscillator, we arrive at the 
effective two-qubit Hamiltonian, 

Heff = ~ A-Ej - - a x ia X 2- (22) 

i 4cos<^o 

The qubit-qubit coupling term in Eq. (|22jl has a clear physical meaning: it results from 
interacting persistent currents in the qubit loops. Indeed, the persistent currents are 
given in terms of Eqs. ©-© by the relation Li = (2e/K)Ejsin(<f>ii), or identically, 

h = -z- Ejsm— -cos0_i = — — — . (23) 

h 2 h d(p +t i 

In the truncated form, this relation reduces to 

J^l^sin^^ (24) 



(neglecting phase fluctuations over the coupling JJ), and the coupling term in Eq. 
can be expressed as an inductive coupling energy of the two persistent currents, L b IiI 2 
with the Josephson inductance of the tilted coupling junction, 

h 2 



4e 2 Ej cos ipo ' 



(25) 
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playing the role of mutual inductance. In the absence of the bias current when the 
JJ potential is not tilted, sin<^ = 0, the persistent currents are not excited and the 
coupling is switched off. When the bias current is applied, the coupling is switched on, 
and its strength increases with the bias current because of increasing persistent currents, 
and also because of decreasing JJ inductance. 

2. 3. Effect of qubit asymmetry 

Let us consider the effect of small a-terms in Eq. ()14p . Although small, the last term 
in this equation proportional to o~ y leads to an interesting qualitative effect in the qubit 
coupling, changing its symmetry. A similar effect is produced by asymmetry of the qubit 
junctions. Although an ideal qubit should consist of identical Josephson junctions, in 
practice the junction parameters may vary at least within the range of a few percent. 
In the asymmetric case, the property of the symmetric qubit to have zero persistent 
current when the bias is turned off is lost. Now a persistent current is spontaneously 
generated, the direction of which depends on the charge state of the SCB. This affects 
the symmetry of the controllable qubit coupling. 

The most important is the variation of the Josephson energy. For an asymmetric 
qubit, the Josephson term in Lagrangian, Eq. (JHJ), has the form, Ej\ cos (fin + Ej 2 cos02i- 
For small junction asymmetry, SEj = Ej\ — Ej 2 <C Ej, the Josephson term in the qubit 
Hamiltonian, Eq. ([lljh acquires the form 

ip 5Ej . if 

- Ej cos— a xi ± ——sm-a yi . (26) 

The second term in this equation, resulting from the junction asymmetry, has the same 
^/-symmetry as the last term in Eq. (jl4j) . They can therefore be considered on the same 
footing and added to the interaction Hamiltonian (|17|) . which now takes the form, 



H int = Eji 



l y?o / y?o SEj 

- sin —{cr x i + cr x2 ) + cos — 2a — a<^ tan — 



[a y2 - o y x / 



.(27) 



.2 2 v ' 2 \ r 2 AEj / 

The additional terms give rise to a small direct qubit coupling of the xy-type, in addition 
to the controllable xx-coupling in Eq. (J22|) . 

^AEjtan^o ^2a - a^otany - ^r^J (a y ia x2 - o- xl a y2 ). (28) 

Although small, this additional coupling term does not commute with the qubit 
Hamiltonian even at the degeneracy point, which may complicate the gate operation 
discussed towards the end of this paper. 



2.4- Residual qubit coupling 

Even in the absence of bias current, and in the symmetric qubit s, there exist small 
circulating currents in the qubit loops because of ground state phase fluctuation in the 
coupling junction. These currents interact via the coupling junction, creating a small 
parasitic coupling of the qubits. The effect is described by the higher order terms 
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neglected in the previous discussion. One term, which does not vanish at (p = is due 
to the interaction via the gate capacitance (cf. Eq. (|27j)). 

H$ f = 2aEj cos(^ /2) 7 K2 - <r y i). (29) 

This term is linear in 7, and it creates a direct parasitic qubit-qubit coupling via the 
mechanism discussed in the previous sections, 

fl2 = WAB,5S!M^. (30) 
cosy9 

This coupling is smaller than the controllable coupling by a factor a 2 = (Cg/iC^) 2 <C 1. 

Obviously, the effect of the junction asymmetry also contributes to this kind of 
residual coupling, and can be included in Eq. (JHDJ) , by making a change, a — > 
a - 5Ej/2Ej. 

Another parasitic term is quadratic in 7, 

^ 2 ] = ^cos^7 2 Ki + M. (31) 

The effect of this interaction is to change the frequency, and hence the ground state 
energy of oscillator (fTUjl. depending on the qubit state configuration. This squeezing 
effect creates a direct residual qubit coupling in the lowest order approximation, 

rr(2) 1 X p fab COS 2 

" -T28 AEj ^(cos^ ) 3 / 2 Xl x2 ' (32) 
This coupling is smaller than the controllable interaction in Eq. (|2*2*|) by a factor, 
hu b /E b j < 1. 



£.5. Maximum coupling strength 

Because of the limitation on the gate operation time imposed by decoherence, it is 
desirable that the qubit coupling is as strong as possible. In our case, the coupling 
strength is generally determined by the parameter A; the strength however increases 
with the applied current bias. This is reflected by a cosine-factor in the denominator 
in Eq. (|22j). which formally turns to zero at (po = tt/2. This corresponds to the point 
when the bias current approaches the critical current value for the coupling JJ. At 
this point the minimum in the tilted Josephson potential disappears, and the junction 
switches to the resistive state, sweeping the qubit phase and thus destroying the qubit. 
Therefore the ultimate limitation on the coupling strength is imposed by the switching 
of the coupling JJ. The latter may even occur at smaller applied current because of 
tunneling through the Josephson potential barrier (macroscopic quantum tunnelling, 
MQT). The assumption of a small MQT rate imposes an additional limitation on 
the coupling strength to the one imposed by the constraints (j2"D|h Indeed, because 
the potential wells of the tilted Josephson potential become shallow with decreasing 
Josephson energy, E 1 } cos ipo, the constraints must be reconsidered, 

Ej < co> bv /cosv?o < -Ejcos^o, (33) 
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clearly putting limitations on the maximum allowed tilt. 

In order to very roughly estimate an upper bound for the maximum coupling 
strength, let us soften requirements (JHBJ), and consider the relations 

Ej ~ uj by J cos tfQ ~ Ej cos if®. (34) 

Both the relations can be fulfilled by applying sufficiently large bias current, and by 
choosing appropriate plasma frequency. The latter can be adjusted by shunting the 
coupling JJ with a large capacitance. The corresponding relations read, 

cos^o~§-, u b ~ \j EjEj . (35) 

The coupling strength for a tilted JJ is given by the phase- dependent coupling parameter 
in Eq. (ED, 

sin 2 (y? /2) 

Kvo) = x — ; • 36 

4 cos <p 

The maximum value of this parameter is estimated by using Eq. (|35|) . 

max \((po) ~ 1 • (37) 

Let us compare this result with the limitation imposed by MQT. For large applied bias 
current, the potential well can be approximated with a cubic curve, and the MQT rate 
is estimated bv|3Ti] 



r /30s _ s ME* (cosyo) 5 / 2 

T M QT = u b \ cos^ e , s=— — . (38) 

V 7r but sin (p 

Suppose that the value u) h ~ (l/2)_E , j(cos<^ ) 5 ^ 2 gives satisfactory small MQT rate 
(~ 10~ 4 Ej according to the following estimates). Under this condition, which is more 
restrictive that the right one in Eq. ([34)1. the relations in Eq. ()35|) become modified, 

cos -(g) 173 , .^ET{E^'\ (39) 

leading to a somewhat smaller maximum coupling parameter, 

(E \ 2/3 

maxA(v? )~ (^J < 1- (40) 
2.6. Charge-phase regime 

So far, we have assumed the qubit island to be in the charge regime Ec 3> Ej, 
where the two lowest charge eigenstates, |n = 0) and |n = 1), serve as the qubit 
basis. However, from an experimental point of view it may be more appealing to 
work in the charge-phase regime Ec ~ Ej because the qubit becomes more stable 
against charge noise when the energy bands flatten [2]. In this regime, the qubit 
states are given by Bloch wave functions, consisting of superpositions of many charge 
states. Nevertheless, as easily seen, the controllable qubit coupling via a current biased 
large JJ will persist also in the charge-phase regime. Indeed, an essential physical 
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characteristic of the qubit- JJ interaction is the persistent current in the qubit loop, Eq. 
([2.2)1 . / = (2e/h) sin(0 + /2) cos</>_. The magnitude of this current is controlled by the 
tilt of the JJ (sin(0 + /2) = sin(</?o/2)), and it is zero when the JJ is idle, regardless of 
whether the qubit is in the charge or charge-phase regime. 

Furthermore, an important property of the charge regime is that the qubit- JJ 
interaction is diagonal in the qubit eigenbasis when the qubit is parked at the charge 
degeneracy point, n g = 1/2, i.e. it has zz-symmetry in this eigenbasis. This property 
simplifies the 2-qubit gate operations discussed later, and it also allows a quantum non- 
demolishing measurement of the qubit by means of current detection using the large JJ, 
as discussed in the next section. We show in this section that this property persists in 
the charge-phase regime. Namely, we show that the SCB Hamiltonian truncated to a 
pair of lowest Bloch states commutes with the truncated current operator at n g = 1/2. 
This leads to direct qubit-qubit coupling of zz-type in the qubit eigenbasis. 

Let us consider the SCB Hamiltonian Hscb, Eq. ©, in the charge basis, \n), and 
separate the part which does not depend on the gate charge, 



E c n(n - l)\n){n\ - Ej (\n + l){n\ + \n - l){n\) 



(41) 



from a small part proportional to the departure from the charge degeneracy point (e.g. 
during single-qubit manipulation) 5n g (t) = 1/2 — n g {t), 



H 2 = 2E c 5n g (t)n\n)( 



n\ 



Hscb — Hi + H 2 . 



(42) 



The notation Ej = 2Ej cos(y?/2) is introduced here for brevity. We split the complete 
set of the charge eigenstates, — oo < n < oo, in the positive and negative charge subsets 
labelled with a =f , |, and m, 1 < m < oo, such that 



m = n, n > 0, 

m = 1 — n, n < 0. 

In the basis |m, a), m = ... ,2,1, the Hamiltonian Hi acquires the form, 
Hi 



H 


-EjU~ 


-EjU 


H \ 



where Hq is tridiagonal, and U contains only a single element, 



Hn 



QE C ~Ej 

-Ej 2E c 

-Ej 

A Hadamard rotation, H, in cr-space, 
1 
72 



-Ej 





u 




1 



H 



(43) 



(44) 



(45) 



(46) 
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takes the basis \m ]), \m |) to |m±) 
matrix in Eq 



Hi = H l- EjUa x ^H 1- EjUa z . 
Then this matrix acquires a block-diagonal form, 



(47) 







#o - EjU 





X 







H + EjU 




X 



(4? 



The above block-diagonal form is suitable for identifying the qubit states in the 
charge-phase regime. Indeed, when the Josephson energy is tuned to zero, the lower- 
corner elements of the blocks, marked with x, correspond to the lowest energy states 
of the SCB. In fact, the Hadamard transformation f)4fij) corresponds to the rotation 
to the eigenbasis of a charge qubit at the charge degeneracy point, which coincides 
with the current eigenbasis. When the Josephson energy increases, the eigenstates of 
the matrix (}4*5|) become superpositions of many charge states, which however does not 
mix the charge superpositions denoted with indices + and — , and can be obtained by 
independent rotations of the matrix blocks. During these rotations, although the two 
lower-corner eigenvalues marked with x do change, they however remain the lowest 
energy levels. This follows from the fact that the eigenvalues of the Mathieu equation 
do not cross when the amplitude of the potential increases [28J. Therefore the charge 
qubit eigenstates develop to the lowest energy Bloch states, which are identified as the 
charge-phase qubit eigenbasis \E + ) and \E_). 

Let us evaluate the form of the current operator, Eq. (|2.2|) . in the charge-phase 
qubit eigenbasis. The current operator in the charge representation is proportional to 
the operator X = \n+ l)(n\ + \n 



X = X Q 1 + Ua x 



l)(n|. In the |mcx)-basis, this operator is written as 

(49) 



where, 



1 1 
1 



(50) 



In the eigenbasis of Hi, i.e. after the Hadamard transformation, this operator acquires 
a block-diagonal form, 



' X + U 








x -u 



X 



which means that X does not couple the states \E + ) and \EJ). 



(51) 
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One implication of this result is that the qubit-qubit coupling in the charge-phase 
regime will still be of zz-type in the truncated Hilbert space, i.e. diagonal in the qubit 
eigenbasis. It also means that current measurement will not mix the qubit states, i.e., 
current detection provides a means for quantum non-demolition measurements. 

Finally we analyze the term H 2 in Eq. (|42|) . which is non-zero only when the gate 



charge deviates from the degeneracy point, 
form, 



H 2 = 2E c 6n g {t) 














In the | racr) representation this term has 



(52) 



where and are diagonal matrices, 

£> T = diag(..., 3,2,1) 
D i = diag(...,-2,-l,0). 

After the Hadamard rotation it acquires the form, 
H 2 -> E c Sn g (t) 



(53) 



(54) 



Thus after truncation to the qubit basis, — provide off-diagonal elements, which 
couple the qubit states and can be employed for qubit manipulation. 



1 









3. Coupling via read-out junctions 

The readout circuit is an important ingredient of the qubit network, which must be 
explicitly included in the consideration. We discuss here the readout method successfully 
tested on a single qubit by the Saclay group [2 a . With this method, the persistent current 
flowing in the qubit loop is excited and measured by using a large Josephson junction 
in the qubit loop, as shown in Fig. |21 To do the measurement, a large dc current is 
applied to the junction so that the net current through the junction either exceeds the 
critical value or not depending on the direction of the persistent current in the loop. In 
the former case, the measurement JJ switches to a resistive state, which is detected by 
measuring a dc voltage across the JJ; in the latter case, no voltage is detected. This 
method of threshold detection is quite invasive, sweeping the phase across the qubit and 
creating a large number of quasiparticles. A recently tested more gentle method [22] 
utilizes an ac driving current with comparatively small amplitude applied to the JJ, and 
measures the qubit-state dependent ac response. 

We now analyze the compatibility of such measurement methods, via a large JJ, 
with our coupling scheme. Before proceeding with the calculations we note that one 
may distinguish two cases: measurement and coupling. In the measurement case, the 
bias current is applied only to one single measurement junction. This will excite the 
persistent current in the corresponding qubit loop, allowing qubit readout, while the 
neighboring qubit loop will remain, as we will see, in the idle state, neglecting the effect 
of a small parasitic coupling, and this qubit will not be destroyed. In the coupling case, 
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bias current sent through both measurement junctions in Fig. |2]will create persistent 
currents in both qubit loops, resulting in the qubit coupling discussed above. This 
physical picture implies that qubit-qubit coupling can be achieved even without sending 
current through the coupling junction 

3.1. Measurement of individual qubits 



Qubit 1 









nl 


Qubit 2 



~r c g 
jv gl 



j v g2 



Figure 2. Applied currents I m i across large measurement JJs can be used for reading 
out the current state of qubit i and for qubit coupling. 



When a measurement J J is included in each qubit loop (Fig. |2J, the corresponding 
terms must be added to the circuit Lagrangian Eq. (j2J): 

h 2 C ■ h 

iu ^mt . o nm t /rr\ 

= 2{2e) 2 ^ mi Ji C ° S 2e miVmi ' ^ ' 

Here (p m i denotes the phase across the measurement junction of the z-th qubit, and J m j 
is the applied current. The phase quantization relations (jSJ) will now change, 

0+,l + V - Vml = 0, 0+,2 - f ~ <Pm2 = , (56) 

giving rise to interaction of the qubits with the measurement JJ, in addition to the 
coupling JJ, in the charge sector as well as in the current sector. As before, it is 
possible, to eliminate the capacitive interaction, but now it is more convenient to do 
it on the Lagrangian level. The interaction via the gate capacitance is eliminated via 
transformation of the qubit variable, 

4>-,i -> 4>~,i + 2a(ip mi =F v 5 ), (57) 

which is equivalent to the transformation in Eq. ()13p (the upper (lower) sign corresponds 
to the first (second) qubit). As already described, this interaction leads to a small 
residual direct yy qubit coupling, which we will omit from the further discussion. 
Similarly, the capacitive interaction via the qubit capacitance C can be eliminated by 
transformation of the measurement JJ variable, 

C 

(pmi -> <Pmi±(3<P, (3 = — — -■ (58) 

^mi ~r Li 

This transformation will only slightly affect the inductive interaction since (3 is small. 



Superconducting qubit network with controllable nearest neighbor coupling 



15 



At this point, we proceed to the quantum description of the circuit and truncate 
the qubit Hamiltonian, assuming the charge regime, Eq 3> Ej. The circuit Hamiltonian 
will take the form, 



where 



H = J2(H i + H mi ) + H JJj (59) 



Hi = — (1 - 2n 9i ) cr zi - Ej cos a xi , (60) 



refers to the qubits, while 

H mi = E™n 2 mi - ETi cos(^ mi ± /3<p) - ^Imi&mi ± P<p), (61) 

and 

Hjj = E b c n 2 -E b jCos<p, (62) 

are the Hamiltonians of the measurement JJ and the coupling JJ, respectively. In these 
equations, E% = {2e) 2 /{2C mi + C), and E b c = {2e) 2 /2C b ^, C 6 e = C b + E l [C{C mi + 
C)/(2C m i + C)]; the current applied to the coupling junction is absent because we will 
focus on the effect of the measurement junctions. 

In the measurement regime only a single external current, say J m i, is applied. The 
steady state point for the 3-JJ network, (<p Q , y^mi.o), is found from Eqs. (J60"|) - (}62|) in the 
main approximation with respect to the small parameters (3 ~ C/C b and A = Ej/Ej, 

h I m \ 



sin (p mlfi 



2e ER 



^° = 2 Sm_ 2~ axl) 

<Pm2,o = Pvo - ^ sin y * 2 " ( 63 ) 

It follows from these equations that indeed the phases across the coupling JJ and the 
second measurement JJ remain negligibly small even though the first measurement 
junction may be biased at the critical level. Thus, the constraint (J2*U|) is essential for 
not disturbing the other qubit while the first qubit is measured. 

3.2. Qubit coupling via readout junctions 

In the case of qubit-qubit coupling, both measurement junctions are biased, while the 
coupling junction is not tilted by external bias, 

siny? mii0 = — — , y? = 0. (64) 

Expanding the potential terms in Eqs. (|5D* j) - (|52*j) around the steady state point up to 
second order with respect to small phase fluctuations, 7 and 7 mi = (p mi — v?mi,cb an d 
neglecting /^-corrections, we may present the Hamiltonian on the form 

2 

H = Y J H l + H osc + H mt . (65) 
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Here 

Hi = ^ (1 - 2n gi ) a zi - Ej cos cr xi , (66) 

is the qubit Hamiltonian, which differs from the one in Eq. by the phase of the 

measurement JJ substituting for the phase of the coupling JJ. The next term, 

H osc = E b c n 2 + E™n 2 ml + E™n 2 m2 + ^(E^ 2 + E™ cos <p ml>ol 2 ml + EJ 2 cos <p m2 , ^ 2 ) (67) 

describes uncoupled linear oscillators, while the interaction is described by the last term, 
XE b 

H int = ~y~ [( A i + A 2) 7 2 + A l7 ^ + A 2l 2 m2 ] + XE b j (A 27m2 - A l7ml ) 7 

+\E b j \{B X - £? 2 ) 7 - B llml - B 2lm2 ] . (68) 

Here we introduced for brevity the following notations, 
Ai = (1/4) cos((p mifl /2)a xi , 

B t = -(1/2) sm((p mij0 /2)a xi . (69) 

Now our goal will be to eliminate the linear terms in gammas in Eq. (J68|) . which can 
be easily done by using oscillator normal modes. To this end we rewrite the potential 
part of Eqs. (ffiTj) . (pj%|) in a symbolic form in terms of a 3- vector 7 = (7, j ml , 7 m2 ), 

±E b jftD + \&)*r + \E b jBr (70) 

Here D is a diagonal matrix representing the free oscillator potentials in Eq. (JBTj) , while 
the 3x3 matrix A and the 3- vector B represent the interaction in Eqs. (J68|) and (|69|) . 
Without loss of generality we may assume the charging energies of the oscillators to be 
equal [3*2] . Then performing rotation to the eigenbasis j' of the matrix D + A A and 
then shifting the variable, 7 = 7' + \D'~ l B' (here the prime indicates a new basis), we 
get, 

l -E b 7D'7 - ^E b B{D + AA)-^. (71) 

The last term in this equation, which is conveniently written in the original basis, gives 
a direct controllable qubit-qubit coupling similar to the one in Eq. fl22j) . 

H int = -XEj sin sin a xl a x2 . (72) 

As expected, this coupling is switched off when one or both measurement junctions are 
idle, and it is switched on only when both the measurement junctions are biased. We 
emphasize that this coupling does not require biasing of the coupling junction. 

The first term in Eq. ()71|) gives, after averaging over the oscillator ground state, the 

oscillator ground state energy, (/^/2)Tr\//5' (remember that D' is diagonal). Treating 
AA as a small perturbation, we find the first perturbative correction to the matrix 
spectrum, D' = D + A diagA. It is easy to see that only the contribution of the 
coupling JJ contains the dependence on the qubit state configuration. The relevant 
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matrix element has the explicit form (l/2)Ej [1 + A(Aj + A 2 )], and yields the residual 
interaction 

1 /iCo> b <£>ml,0 ^m2,0 / 7Q \ 

= - ^ * J ~Ef cos ~~ ^T - cos ~~ ^T - o"xi^2- (73) 
This is a small residual interaction substituting for Eq. (J3*2|) in the present case. 



4. Multiqubit network 

To implement useful quantum algorithms, controllable systems with large numbers of 
qubits are needed. In this section we will show that the effective qubit-qubit coupling 
derived in sections 121 and H3 can be generalized to a chain of iV qubits with each qubit 
being coupled to its nearest neighbors via current-biased JJs and each having its own 
read-out device, as shown on Fig. El 




Figure 3. A system of N coupled charge qubits. V g i controls individual qubit, whereas 
Ibi controls the coupling of qubits i and i + 1. 7 m , can be used to read out the current 
state of qubit i and also for qubit coupling. There are no coupling JJs at the ends, 

lfi = ifiN = 0. 



4-1. Circuit Hamiltonian 

The Lagrangian of the iV-qubit circuit presented in Fig. 01 can be written as a 
straightforward generalization of Eqs. (J2j), 0, (jlj), and (J55f . 

N N-l 

£ = £ [L S CB,i + L mi ] + £ L JJ:i . (74) 

i=l i=l 

Since now there are two coupling JJs per qubit loop, the flux quantization relation (j3Sj) 
must be extended, 

+ (fi- 4> mi - (fii-l = 0, (75) 
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leading to a more complex form of the interaction among the qubits and the coupling and 
measurement JJs. Nevertheless, our previous strategy for elimination of the interaction 
in the charge sector still works. Generalizing Eq. (|57j) . 

4>-,i -> 4>-,i + 2a (Vmi + <-Pi-l + Vi)- (76) 

allows us to decouple qubit charges; the resulting weak interaction in the current sector 
yields a direct parasitic yy qubit coupling, similar to the one in Eq. (fHUJ) . Further 
transformation, generalizing Eq. (|58j). 

C 

fmi -> Vmi + P{<fi - (pi-l), P = —77, (77) 

decouples the charges of the measurement JJs, and yields weak additional interaction 
in the current sector, which will be also omitted. 

After the transformations (ffT)|) and (|77|). the only capacitive interaction which 
remains in the Lagrangian (J73j) is the interaction among the coupling JJs. This 
interaction is proportional to a small qubit capacitance C, while the diagonal terms 
are proportional to much larger capacitances of the coupling JJs, Cm ^> C. We assume 
that the J J capacitances are different so that C^ — C7jy ^> C; this is a realistic assumption 
because a spread of the junction characteristics during fabrication usually exceeds 10%. 
Under this assumption, the diagonalization of the capacitance matrix will introduce 
small corrections to the inductive interaction, corrections which will not provide any 
qualitative changes and can be omitted. 

With diagonal kinetic terms in the Lagrangian (|74p. it is straightforward to proceed 
to the truncation of the quantum Hamiltonian, 

N N-l 

H = J2[H l + H mi \ + ]T H JJti . (78) 
i=i i=i 

In the charge regime, Ej <C Ec, the qubit Hamiltonian will take the form, 

Hi = — (1 - 2n gi )a zi - Ej cos a xi , (79) 

while the large- JJ terms will not change, 

m 2 m h 

H m i = E Ci n mi — Ej i COS<£> m j — —Imi^mii 

H JJ:i = E h c n\ - E b j cos ^ - ^hiVi- (80) 

In these equations, the effective capacitances of the measurement JJs are the same 
as in Eq. (|65jh while the effective capacitances of the coupling JJs are straightforward 
generalizations of the one in Eq. (JfSJ), namely = Cb + C[{C mi + C)/{2C mi + C) + 
(Cm,i+1 + C)/(2C mj i + i + C)j. 

Direct qubit- qubit coupling 



The next step in the derivation of the direct qubit-qubit coupling is to eliminate the 
large JJs, following the previous procedure for the two-qubit case. After expanding the 
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Hamiltonian, Eqs.(|78 |) -(|8U )) . with respect to small fluctuations around the steady state 
points, ((fio, Prnifi), determined by the applied controlling and measurement currents, 

hi 
2e~E~ b , 



sm (p mifl 



2eE 



smcfio = —-^ , 



11) 



Ji 



we get the qubit terms (ffD^l with steady state phases, (JHTj) . in the Josephson terms. The 
qubits interact with a subnetwork of the linear oscillators, 



N-l 

E 

i=l 



E 



N 



+ E 

i=l 



E 

E Ci n mi + Y"COS0 miO T. 



,(82) 



via the interaction Hamiltonian, which also connects the oscillators, 



H, 



XEj 7V ~ 1 



int 



E [(A t + A m ) 7i 2 -2A 



i=l 



iTiTi— 1 



- JV JV-1 

E A ^L + 2E( A 

.i=l i=l 
AT-1 iV 

^ (S< - J B i+1 )7 i - Si7„ 
1=1 1=1 



53) 



The quantities A and B now contain also the phases of the two coupling JJs (cf. Eq. 



A 1 PiO ~ 4>mi0 ~ <Pi-l,0 

Aj = - cos a xi , 

Bi = - sm a xi . 



J4) 



The interaction (|83|) can be presented in the symbolic form of Eq. (J7Uj) by 
introducing the 2N — 1-vector 7 = (71, 72, ... , 7mi, 7m2, ■ • the 2N — 1 x 2N — 1 matrix 
A representing the oscillator interaction, and the 2N— 1-vector B representing the qubit- 
oscillator interaction. Then we proceed to Eq. ([71 p by performing the diagonalization, 
and shifting the oscillator variables as described after Eq. (J7UJ). The result of this 
procedure is as follows [22] : Assuming no applied measurement currents, the coupling 
induced by only tilting coupling JJs has the form, 



H, 



int 



N-l 

E 



A-Ej . <PiO — <Pi-X,Q ■ Vi+1,0 — PiO 

sm - sm '— cr xi a x i+ i. 



t r{ 4cosy2 i0 2 2 

On the other hand, when the coupling JJs are kept idle while the measurement junctions 
are biased, the coupling has the form, 



H, 



N ^XE 



int 



E 

i=l 



J . (PmiO ■ (PmA+1,0 

- sm — — sm o- xi a X A +1 , 



(86) 



Whichever way the coupling is initiated, one is allowed to simultaneously perform a 
number of two-qubit gates on different qubit pairs, as long as the qubit pairs are 
separated by at least one idle qubit. The small residual xx coupling resulting from 
the shift of the oscillators ground energy is restricted to the neighboring qubits and 
given by Eq. 
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We conclude this section with a discussion of the effect of different Josephson 
energies of the qubit junctions, E Jiy and the coupling JJ, E b M . This variation can 
easily be taken into account by introducing numerical scaling factors, Eji = £iEj, and 
E h M = i\E b j. Then, while deriving Eq. ([53]). these scaling factors can be included in the 
definition of the quantities A, and Bi in Eq. As a result, the coupling energies 

XEj in the final results, Eqs. (JSSJ), and (JBSJ), are replaced by 



\Ej = § - (87) 
E b j E b M 1 ; 



5. Gate operations with the qubit network 

All quantum algorithms can be implemented using a limited universal set of gates. 
One such set consists of the controlled-NOT (CNOT) gate together with single qubit 
gates In this section, we will describe how to perform a CNOT gate on two 

neighbouring qubits in the above mentioned charge qubit network. Using a sequence 
of two-qubit operations on nearest neighbours only, two-qubit operations on arbitrary 
qubits in the chain can be performed ^2] • The CNOT gate presented here is composed 
of a control-phase (CPHASE) gate and two kinds of single qubit gates, a phase gate 
and the Hadamard gate. 

By default, during qubit operations the qubits are parked at the charge degeneracy 
point, where they are more stable against charge noise and the qubit levels are 
maximally separated from higher states. The computational basis is chosen to be the 
current basis, which is the eigenbasis at the charge degeneracy point and differs from 
the charge basis by the rotation a x <-> o z . 

5.1. Single qubit operations 

When I bi = h,i-i = l m % = 0, qubit i is disconnected from the network to first order and 
single qubit gates can be performed. The time evolution is determined by the single 
qubit Hamiltonian H iy Eq. ([To*)) or (JMJ), 

Hi = E c Sn gi (t)a xi - Eja zi , (88) 

here 5n g i(t) = 1/2 — n gi (t) is the deviation from the charge degeneracy point. 

In the idle state, the non-zero energy level splitting results in a phase gate being 
performed on the qubit; 

_ J |0) - e^\Q) 
he -\ |1> - e-^\l) ' (89) 

where 9 depends on the elapsed time T through T = 6/{2Ej). 

A particularly useful case is the S , 37r / 2 -gate which will be referred to as the Z-gate, 
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(91) 



Another useful single qubit operation is the Hadamard gate H, 
H f |0> - ^(|0> + |1» 

1 |i> - ^(|o)-|i» ' 

which can be implemented by applying a microwave pulse at the gate [HUI2|, Sn gi (t) = 
Acos(2Ejt), during a time T = ir/2A. Choosing the amplitude A involves a trade-off 
between keeping the operation time short and minimizing the deviations from the charge 
degeneracy point. 



5.2. Two-qubit gates 

A two-qubit gate involving qubits i and i + 1 is created by applying a bias current 7^ 
at the intersection between the two qubits, or by simultaneously applying measurement 
currents I m i and I m ,i+i- The qubits are coupled according to the coupling terms Eqs. 
(I85j) and (|86jl . while their individual time evolutions are determined by Hi, Eq. (fTH|) or 
f!66|) . As an example, when applying the bias current the Hamiltonian of the two 
interacting qubits reads, 



Hi + H i+ i + H, t 



(i) 

int 



-Ej cos— [a zi + cr z ,i+i) 



■ 2 V'iO 

sin —a zi a Z) i +1 . 

4 cos (fio 2 



(92) 



CPHASE 



(93) 



Choosing operation time and bias current amplitude properly results in the entangling 
control-phase (CPHASE) gate, 

|11) -> i|ll) 
|10) - |10) 
|01) -> |01) 
{ |00) ^ z|00) 

Moreover, a CNOT gate is created by combining CPHASE with single qubit gates such 
as the Z-gate, Eq. (|9T)|) . and the Hadamard gate, Eq. (jHH), as shown in Fig. HJ Thus 
it is possible to perform a universal set of quantum gates, and therefore any quantum 
algorithm, with the investigated charge-qubit network. 




CPHASE 



QUBIT 1 



QUBIT 2 


H 






C 





z 




H 




) 







Figure 4. A CNOT operation using single qubit gates (Hadamard and phase gates) 
and CPHASE. Time runs from left to right. 
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